Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Variable gamma-ray emission from the Be/X-ray transient 

A0535+26? 



o 
o 



O 
(N 



O 



43 

Or 

6 



% 



Gustavo E. Romero*, M.M. Kaufman Bernado, Jorge A. Combi*, Diego F. Torres 

Institute Argentino de Radioastronomfa, C.C.5, (1894) Villa Elisa, Buenos Aires, Argentina 

Received / Accepted 

Abstract. We present a study of the unidentified gamma-ray source 3EG J0542+2610. This source is spatially 
superposed to the supernova remnant G180.0-1.7, but its time variability makes unlikely a physical link. We have 
searched into the EGRET location error box for compact radio sources that could be the low energy counterpart 
of the gamma-ray source. Although 29 point-like radio sources were detected and measured, none of them is 
strong enough as to be considered the counterpart of a background gamma-ray emitting AGN. We suggest 
that the only object within the 95% error box capable of producing the required gamma-ray flux is the X-ray 
transient A0535+26. We show that this Be/accreting pulsar can produce variable hadronic gamma-ray emission 
through the mechanism originally proposed by Cheng & Ruderman (1989), where a proton beam accelerated in 
a magnetospheric electrostatic gap impacts the transient accretion disk. 
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1. Introduction 

The EGRET instrument of the late Compton Gamma Ray 
Observatory detected 271 point-like gamma-ray sources at 
energies E > 100 MeV during its lifetime (Hartman et al. 
1999). The majority of these sources still remain unidenti- 
fied. Excluding 6 possible artifacts, there are 75 sources at 
low galactic latitudes \b\ < 10° still not conclusively iden- 
tified with objects seen at lower energies (Romero et al. 
1999). Several of these low-latitude, presumably galactic 
sources seem to display significant levels of variability over 
timescales of months (Tompkins 1999, Torres et al. 2001a, 
b). A few particular cases of variable gamma-ray sources 
at low latitudes have been recently studied by Tavani et 
al. (1997, 1998), Paredes et al. (2000) and Punsly et al. 
(2000), but there is not yet consensus on their nature. 
The variety of behaviours displayed by the sources seems 
to suggest that there exist more than a single class of 
gamma-ray emitting objects in the Galaxy. To establish 
the nature of these peculiar objects is one of the most 
urgent problems of high-energy astrophysics. 

In this paper we present a study of the gamma-ray 
source 3EG J0542+2610. We shall show that the only 
known object within the 95% confidence location contour 
of the source capable of generating the observed gamma- 
ray emission is the Be/X-ray transient A0535+26. We 
shall argue that the gamma-rays are produced during the 
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accretion disk formation and loss phases in each orbit, 
through hadronic interactions between relativistic protons 
accelerated in an electrostatic gap in the pulsar magneto- 
sphere and the matter in the disk. The structure of the 
paper is as follows. In the next section we summarize the 
main characteristics of 3EG J0542+2610 and the results 
of our analysis of the radio field within the EGRET loca- 
tion error contours. We then present the peculiarities of 
the system A0535+26. In Section 4 we present a model for 
the gamma-ray production based on the accelerator gap 
model by Cheng & Ruderman (1989, 1991). The main dif- 
ference between our treatment and the accreting pulsar 
gamma-ray model of Cheng et al. (1991) lies in the in- 
troduction of time variability through the variation of the 
matter content of the disk, which is a transient structure 
according to recent observations. 

2. The gamma-ray source 3EG J0542+2610 and 
the surrounding radio field 

The best estimated position of the gamma-ray source 3EG 
J0542+2610 is at (I, 6)«(182.02 -1.99). Its 95% confi- 
dence location contour overlaps with the shell-type su- 
pernova remnant (SNR) G180.0-1.7 (Romero et al. 1999, 
Torres et al. 2001c). The gamma-ray flux for the combined 
EGRET viewing periods is (14.7±3.2) 10" 8 ph cm~ 2 s -1 . 
Since the flux is highly variable on timescales of months, 
it should be originated in a compact object and not in 
the extended SNR. In Figure 1 we show the EGRET light 
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curve. The source switches between periods of clear de- 
tections and periods when only upper bounds to the flux 
can be determined. The variability analysis of Torres et 
al. (2001b) assigns to 3EG J0542+2610 a variability in- 
dex / = 3.16, which means that the source variability 
level is 4.32er above the average (spurious) variability of 
all known gamma-ray pulsars. Tompkins' (1999) variabil- 
ity index for this source (r = 0.7) also indicates that the 
source is variable. 

If a pulsar origin is discarded due to the high vari- 
ability and the steep spectral index (r = — 2.67 ± 0.22), 
we are left with two main possibilities: 1) the source is 
a background, unnoticed gamma-ray blazar seen through 
the galactic plane, or 2) it is a galactic compact object 
with an energy budget high enough as to generate sig- 
nificant gamma-ray emission and, at the same time, has 
not the stable properties usually associated with isolated 
pulsars. 

In the standard model for gamma-ray blazars, the high 
energy photons are produced by ultra-relativistic electrons 
through inverse Compton interactions with seed lower en- 
ergy photons that can come from external sources (e.g 
the accretion disk) or from the jet itself (see Krolik 1999 
and references therein). The model is supported by the 
fact that all identified gamma-ray AGNs are strong flat- 
spectrum radio sources (Mattox et al. 1997). 

In Figure 2, lower panel, we show a 1.4-GHz VLA map 
made with data from the NVSS Sky Survey (Condon et al. 
1998), where all point-like radio sources within the 68% 
confidence contour of 3EG J0542+2610 can be seen. The 
measured characteristics of these 29 sources are listed in 
Table 1. Most of them have no entry in any existing point 
source catalog. In those cases where we were able to find 
positional counterparts at other frequencies we have esti- 
mated the spectral indices, which are also shown in the 
table. Most of these sources are very weak, at the level 
of a few mjy. No strong (at Jy level), flat or nonthermal 
source is within the location error box of the gamma-ray 
source. The strongest radio source (No. 23 in our table) 
has a rather steep spectrum and is a factor ~ 10 below the 
minimum flux density of firm gamma-ray blazar identifi- 
cations given by Mattox et al. (1997). The nature of this 
source is not clear at present; it could be a background 
weak radio quasar. The fact that it is not seen at X-rays 
seems to argue against a galactic microquasar or any other 
kind of accreting source. 

At X-ray energies the most significant source within 
the EGRET error box is the X-ray transient A0535+26, 
which is discussed in the next section. This source does 
not present significant radio emission and consequently it 
cannot be seen in our maps. We have indicated its po- 
sition with a star symbol in Fig. 2, middle panel. We 
also show in this figure the direction of the proper mo- 
tion of the system, as determined by Lee Clark & Dolan 
(1999). The upper panel shows the entire radio field as 
determined from radio observations with Effclsbcrg 100- 
m single dish telescope at 1.408 GHz (data from Reich et 
al. 1997). The middle panel present an enhanced image ob- 



Table 1. Point radio sources within the inner loca- 
tion probability contours of the gamma-ray source 3EG 
J0542+2610 



Number 


I 


b 


Si. 42 GHz 


Spectral 


(on the map) 


(deg) 


(deg) 


(mjy) 


index 


1 


181.46 


-1.56 


18.00 




2 


181.50 


-2.18 


13.15 


„ 0.408 o O 

«i.42 = -2.9 


3 


181.56 


-1.71 


62.90 


alif = -1.3 


4 


181.59 


-1.75 


78.67 


a°df = -1.6 

„ 1.42 n n 
04.85 = -0.9 


5 


181.61 


-2.22 


7.65 




6 


181.62 


-2.09 


7.39 




7 


181.67 


-1.88 


22.47 




8 


181.68 


-2.39 


14.61 




9 


181.76 


-1.88 


25.00 




10 


181.78 


-2.22 


35.00 


„ 1.42 n o 
«4.85 = 0.2 


11 


181.81 


-1.58 


27.00 




12 


181.84 


-1.56 


18.22 




13 


181.95 


-2.35 


11.80 




14 


182.02 


-1.89 


20.32 




15 


182.05 


-1.80 


53.40 




16 


182.18 


-2.01 


37.30 




17 


182.18 


-1.40 


33.00 




18 


182.27 


-2.22 


8.70 




19 


182.27 


-1.73 


18.88 




20 


182.32 


-1.71 


9.91 


a?;!! 1 = -2.27 


21 


182.33 


-2.35 


10.30 




22 


182.35 


-2.38 


6.00 




23 


182.38 


-1.72 


225.40 


agifi = -0.86 
all? = -0.75 


24 


182.40 


-1.83 


18.72 




25 


182.41 


-2.41 


20.67 




26 


182.43 


-2.44 


11.90 




27 


182.44 


-2.53 


22.33 




28 


182.50 


-2.48 


57.80 


alif = -1.33 


29 


182.58 


-1.59 


41.86 





tained at 2.695 GHz with the same telescope (Fuerst et al. 
1990), where the gamma-ray location probability contours 
have been superposed (Hartman et al. 1999). We have 
processed these large-scale images using the background 
filtering techniques described by Combi et al. (1998). 

3. The X-ray transient A0535+26 

A0536+26 is a Be/X-ray transient where the compact ob- 
ject is a 104s pulsar in an eccentric orbit around the B0III 
star HDE 245770 (Giovanelli & Sabau Graziati 1992). Be 
stars are rapidly rotating objects which eject mass irreg- 
ularly forming gaseous disks on their equatorial planes. 
If there is a compact companion in a close orbit, accre- 
tion from the star can result in strong X-ray emission. In 
the case of A0535+26, strong and recurrent X-ray out- 
bursts are observed with a period of 111 days, which has 
been identified with the orbital period (Giovanelli & Sabau 
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Table 2. Physical parameters for A0535+26 (from Janot- 
Pacheco et al. 1987 and Giovanelli & Sabau Graziati 1992) 



01/01/91 01/01/92 01/01/93 01/01/94 01/01/95 01/01/9S 

Viewing Period 

Fig. 1. Flux evolution of 3EG J0542+2610 through single 
viewing periods. 

Fig. 2. Radio maps at different resolutions of the field con- 
taining the gamma-ray source 3EG J0542+2610. Upper 
panel: 1408-MHz image obtained with Effclsberg 100-m 
single dish telescope. The shell-type SNR G180.0-1.7 is 
clearly visible in the map. Contours are shown in steps of 
25 mJy beam -1 , starting from 20 mJy beam -1 . Middle 
panel: 2695-MHz map obtained by the same telescope. 
Contours in steps of 30 mJy beam -1 , starting from 30 
mJy beam -1 . The position of A0535+26 is marked by 
a star symbol. The arrow indicates the direction of the 
proper motion. EGRET confidence location contours are 
superposed to the radio image. Lower panel: VLA im- 
age of the inner region at 1.4 GHz. Contours in steps of 
1 mJy beam -1 , starting from 1 mJy beam -1 . The diffuse 
background emission has been removed from the first two 
maps. 



Graziati 1992). It is generally agreed that these outbursts 
occur when the accretion onto the neutron star increases 
at the periastron passage. The average ratio of the X-ray 
luminosity at the periastron to that of the apoastron is 
~ 100 (Janot-Pacheco et al. 1987). In Table 2 we list the 
main characteristics of the A0535+26 system. 

During a major outburst of A0535+26 in 1994, 
the BATSE instrument of the Compton Gamma Ray 
Observatory detected a broad quasi-periodic oscillation 
(QPO) in the power spectra of the X-ray flux (Finger et 
al. 1996). The QPO component was detected during 33 
days, with a central frequency that was well correlated 
with both the hard X-ray flux and neutron star spin-up 
rate inferred from pulse timing. Finger et al.'s (1996) ob- 
servations are the first clear evidence that an accretion 
disk is formed during giant outbursts. 

Using the simultaneous variations of the spin and the 
QPO frequencies in the context of the beat frequency 
model (Alpar & Shaham 1985), Li (1997) has determined 
the evolution of the ratio £ of the inner accretion disk ra- 
dius to the Alfven radius. He found that during the initial 
rise of the outbursts £ quickly increased from ~ 0.6 to 



Primary spectral type 


B0III 


Primary mass 


9 - 17 M @ 


Secondary mass 


< 2.7 M 


Primary mass loss rate 


7.7 10 -7 M Q yr -1 


Distance 


2.6 ± 0.4 kpc 


Lx peak (average) 


7.5 ± 2.4 10 36 erg s" 1 


X-ray pulse period 


104 s 


Orbital period 


111 ±0.5 d 


Orbital eccentricity 


0.3 - 0.8 


j max / r mm 


~ 100 



~ 1, indicating a transition of the accretion process from 
spherical accretion before the outbursts to disk accretion 
during the high X-ray luminosity phase. 

Prior to the direct evidence for a transient accretion 
disk in A0535+26, Motch et al. (1991) have already sug- 
gested, on the basis of an analysis of the long-term X- 
ray, UV, and optical history of the system, that two dif- 
ferent types of interactions exist between the Be star 
and the pulsar, accordingly with the dynamical state of 
the highly variable circumstellar envelope. "Normal" out- 
bursts would occur for high equatorial wind velocities 
(~ 200 km s -1 ), whereas "giant" outburst would result 
from lower equatorial wind velocities (~ 20 — 80 km s -1 ) 
which allow the formation of a transient accretion disk. 

It is interesting to notice that, during the initial 
stage of the 1994 outburst, when QPOs were found in 
A0535+26, the gamma-ray source 3EG J0542+2610 was 
detected by EGRET with a flux of (39.1 ± 12.5) 10 -8 ph 
cm -2 s -1 (viewing period 321.1: February 8 - 15, 1994). 
But when A0535+26 was at the peak of its X-ray lumi- 
nosity, on Ferbuary 18, the gamma-ray source was not de- 
tected (viewing period 321.5). The gamma-ray emission 
seems to have been quenched precisely when the accre- 
tion disk was well-formed and maximally rotating. In the 
next section we present a model that can account for the 
gamma-ray production in A0535+26 necessary to explain 
the EGRET source 3EG J0542+2610 and that is in agree- 
ment with our present knowledge of the Be/X-ray tran- 
sient source. 



4. Model 

Our purpose in this section is to show that there exist 
a plausible mechanism that could explain the gamma-ray 
emission of the source 3EG J0542+2610 as originated in 
the X-ray transient A0535+26. This mechanism should be 
capable of predicting the observed gamma-ray flux, the 
variability in the lightcurve, the fact that no gamma-ray 
emission was observed on February 18, when A0535+26 
was at the peak of the X-ray outburst, but also that it was 
positively detected in the previous days when the X-ray 
flux was rising, and finally the fact that A0535+26 is not a 
non-thermal radio source. This latter restriction seems to 
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B*0 




Gauss, respectively. The parameter (3 = Ir^jr^ is twice 
the ratio of the inner accretion disk radius to the Alfven 
radius. For the particular case of A0535+26 we can adopt 
(3 ~ 1, according to the estimates by Li (1997). 

Protons entering into the gap are accelerated up to en- 
ergies above E p ~ eAV„i ax whereas inverse Compton and 
curvature losses would limit the energy gain of electrons 
and positrons to lower values. The maximum current that 
can flow through the gap can be determined from the re- 
quirement that the azimuthal magnetic field induced by 
the current cannot exceed that of the initial field B (Cheng 
& Ruderman 1989): 



Fig. 3. Sketch of the magnetosphere model for A0535+26 
when 0<j > ^*- The gray region co-rotates with the star, 
whereas the hatched region co-rotates with the accretion 
disk. An empty plasma gap where E • B ^ exists around 
O* • B = 0. Protons from the stellar surface are acceler- 
ated in this gap and collide with the disk, where they 
produce 7r°-decay gamma-rays. Adapted from Cheng & 
Ruderman (1991). 

suggest a hadronic origin for the gamma rays. Otherwise, 
relativistic electrons should also produce synchrotron ra- 
dio emission. 

Cheng & Ruderman (1989, 1991) have studied the dis- 
turbances produced in the magnetosphere of an accreting 
pulsar when the Keplerian disk rotates more rapidly than 
the star. As in the case of equal angular velocities, when 
fi* < fid, inertial effects of electrons (— ) and ions (+) lead 
to a complete charge separation around the "null surface" 
$1, • B = 0. However, since now the equatorial plasma be- 
tween the inner accretion disk radius tq and the Alfven 
radius rj± co-rotates with the disk, whereas the rest of 
plasma co-rotates with the star, an electrostatic gap with 
no charge at all is created around the "null surface" (see 
Figure 3). In this gap E'B /0 and a strong potential 
drop is established (see Cheng & Ruderman 1991 for de- 
tails). 

Cheng & Ruderman (1989) have shown that the po- 
tential drop along the magnetic field lines through the gap 
is: 



AV m 



B s R 3 ft d {r ) 



r c 



4 io i4 /r b ^ -f- 



M 



1/7 



R, 



-^L^B-r V, (1) 



where B s is the neutron star's surface dipole magnetic 
field, R is its radius, tq is the inner accretion disk radius, 
M is the compact star mass, and L37 is the X-ray luminos- 
ity in units of 10 37 erg s _1 . The radius and magnetic field 
in the second expression are in units of 10 6 cm and 10 12 



J„ 



cB s R 3 r t 
1.5 10 24 /3 





24 o-2 



M 

M© 



-2/7 



pl/7r4/7 R -l/7 1 

Kg -L37 13^ GSU S 



(2) 



As it can be seen from Fig. 3, the proton current flux 
is directed from the polar cap of the star, where the ac- 
creting material impacts producing strong X-ray emission 
and abundant ions to be captured by the gap, towards 
the accretion disk. Electrons move in the opposite sense 
in order to cancel any net charge current flow, keeping in 
this way E.B ^ in the gap. 

The collision of the relativistic proton beam into the 
disk will produce hadronic interactions with copious ir° 
production (Cheng et al. 1990). These it will quickly de- 
cay into 7-rays that could escape only if the disk density 
is sufficiently low. Otherwise, they will be absorbed by 
the matter and then re-emitted as X-rays. The interac- 
tion of relativistic protons with a thin hydrogen layer has 
been studied by Cheng et al. (1990) in the context of their 
model for the Crab pulsar. The 7r°-decay 7-rays will only 
escape insofar as the column density of the disk would 
be £ < 100 g cm" 2 . From Novikov & Thorne (1972), we 
have: 



7 io 2 /r 3/5 



M 



-3/35 



D -30/35 r 27/35 n 
n 6 "^37 a 



12/35 
12 S Cm 



(3) 



In a system like A0535+26, the column density of the 
disk will evolve with time following the variations in the 
X-ray luminosity. All other physical parameters in Eq. (pi) 
remain constant along the orbital period. If at some point 
near the periastron, when the X-ray luminosity is close 
to its maximum, the column density exceeds the value 
~ 100 g cm" 2 , the medium will be no longer transparent 
to gamma-rays and the gamma-emission will be quenched 
(we are assuming that the disk is seen from the opposite 
side to that where the current impacts). This could explain 
the fact that A0535+26 has not been detected by EGRET 
at the peak of its X-ray outbursts on February 18, 1994. 
The detection, instead, was clear a few days before, during 
the previous viewing period. 
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The expected hadronic 7-ray flux from A0535+26 on 
Earth will be: 



F(E 7 ) 



1 



00 ^2 



d 2 N p {E 7T o) 
AtlD 2 J E dt dE v o 



dE p , 



(4) 



where D is the distance to the source, AJ7 is the beaming 
solid angle and: 



d 2 N p (E 7r o) _ 2 J max X < m„o > da(E pi E^a) 



dt dE~o 



dE„ 



(5) 



is the differential production rate of neutral pions from pp 
interactions. In this latter expression, da(E p , E 7V o)/dE 7T o 
is the differential cross-section for the production of ir°- 
mesons of energy E n o by a proton of energy E p in a pp 
collision, and < m^a > is the mean multiplicity of 7r°. 

For calculation purposes we shall adopt the cross sec- 
tion given by Dermer (1986a, b) and a beaming factor 
Aft/Air w 0.3 as used by Cheng et al. (1991). We shall 
evaluate the total gamma-ray luminosity between 100 
MeV and 20 GeV of A0535+2610 at an early epoch t 
of the disk formation, when the X-ray luminosity ratio 
is L^ ax /L^- ~ 10. If the absorption feature observed near 
110 keV is a cyclotron line, the polar magnetic field of 
the neutron star results 9.5 10 12 G (Finger et al. 1996). 
In our calculations we adopt this value along with /3 ~ 1, 
Rq = 1 and M = 1.4 M Q (Janot-Pacheco et al. 1987). 
We then obtain that the column density of the disk is 
~ 42.5 g cm -2 at this stage and, consequently, gamma- 
rays are not absorbed in the disk material. The potential 
drop in the electrostatic gap results AV max ~ 2.5 10 13 V 
whereas the proton current deposited into the disk is 
N p = J m ax/e ~ 4.8 10 32 s^ 1 . Using Eq. (|) for photons in 
EGRET's energy range we obtain: 

F(100 MeV < E 1 < 20 GeV) « 5.1 10~ 8 ph cm" 2 S _1 .(6) 

This flux is consistent with the lower EGRET detec- 
tions. As the accretion onto the neutron star increases the 
gamma-ray flux also increases, reaching a maximum when 
X « 100 g cm" 2 . At this point the flux has grown by a fac- 
tor ~ 4 (notice the dependence on J max in addition to X) 
and photon absorption into the disk becomes important 
quenching the radiation. The gamma-ray source, then, is 
not detected when the peak of the X-ray luminosity oc- 
curs. The additional gamma-ray flux expected from sec- 
ondary electrons and positrons radiating in the the strong 
magnetic fields anchored in the disk can explain even 
higher EGRET fluxes (Cheng et al. 1991). A more detailed 
analysis in this sense will be presented elsewhere, includ- 
ing spectral considerations. The observed spectrum de- 
pends on the proton injection spectrum at the gap. Strong 
shocks near the polar gap can accelerate the protons up 
to a power law which will be preserved through the prop- 
agation across the gap and imprinted into the ir° gamma- 
ray spectrum. Additionally, e + e~ cascades induced in the 
disk material can result in a relativistic bremsstrahlung 
contribution. Electromagnetic shower simulations are in 
progress in order to determine this contribution (Romero 
et al., in preparation). 



5. Discussion 

The model here outlined does not imply strict peri- 
odicity because of the highly chaotic nature of the 
Be stellar winds, which can significantly vary on short 
timescales producing strong changes in the accretion rate. 
Notwithstanding, the general prediction that the peak of 
the gamma-ray emission should not be coincident with 
the maximum X-ray luminosity during a given outburst 
can be used to test the general scenario proposed here. 
At present, the poor time resolution of the gamma-ray 
lightcurve does not allow correlation studies. EGRET data 
for the best sampled X-ray outburst (in February 1994) 
consist only of two viewing periods, as it was mentioned, 
and in one of them the source was not detected. In the fu- 
ture, however, new instruments like GLAST could provide 
the tools for these kind of investigations. 

TeV emission should be produced in the accretion disk 
according to Eq. (ph, although degradation effects during 
the propagation in the strong magnetic and photon fields 
around the accretion disk could suppress much of it. The 
magnetic field at the inner accretion disk radius is (e.g. 
Cheng et al. 1991): 



-3 R -5/7 R -9/7 r 6/7 



B(r ) =3 10* /T d B^"R-""V 37 



(M/M )" 3 / 7 G (7) 



which for the typical values of A0535+26 yields B ~ 10 4 
G. TeV gamma-rays can be absorbed in the field through 
one-photon pair production above the threshold given by 
(e.g. Bednarek 1993): 



X 



E 1 Bsu\6 
2m P c 2 B cr 



1 
15' 



(8) 



where is the angle between the photon momentum and 
the magnetic field and B cr is the critical magnetic field 
given by B CI = m 2 c 3 /eh ~ 4.4 10 13 G. This process, then, 
suppress gamma-ray photons with energies higher than 
~ 300 TeV in A0535+26. Since this value is well above 
the energy of the protons that impact on the disk, we find 
that absorption in the magnetic field should not occur. 

However, TeV photons with lower energies should be 
absorbed by two-photon pair production in the accretion 
disk X-ray photosphere. Calculations by Bednarek (1993) 
show that the optical depth quickly goes to values above 
1 for disks with luminosities L37 ~ 1. The opacity effects 
can reach even GeV energies leading to a steepening in the 
spectrum respect to what is expected from a pure pion- 
decay mechanism. The fact that the observed spectrum in 
3EG J0542+4610 has an index T ~ -2.7 seems to sup- 
port the idea that important absorption is occurring in 
the X-ray photosphere of this source. Additional spectral 
modifications should be produced by the emission of sec- 
ondary pairs in the magnetic field close to the disk (Cheng 
et al. 1991, Romero et al., in preparation). 

In their original model, Cheng & Ruderman (1989) 
suggested that if the disk is sufficiently dense then 
the gamma-rays could be produced only in a moving 
low-density "window" . This window would collimate a 
pencil beam of gamma-rays aligned with the magnetic 



(i 
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axis. Consequently, a pulsed emission could exist with 
the same period of the X-ray source (104 sec in the 
case of A0235+26). Small changes in the gamma-ray 
period might be produced by the radial motion of the 
"window" . It would be interesting to test whether this 
pulses are present in A0235+26. However, since the disk 
is a transient structure and only would excess the critical 
density during a few days per orbit, the number of photon 
counts in EGRET data are too low to allow a periodicity 
analysis as in the case of isolated and stable gamma-ray 
pulsars. Instruments with higher sensitivity like GLAST 
could sum up over several viewing periods in order to 
look for these features. 



6. Summary 

We have shown that the Be/X-ray transient system 
A0535+26 can be also a transient gamma-ray source un- 
der very reasonable assumptions. The existence of a QPO 
phenomenon detected by BATSE during the 1994 X-ray 
outburst provided direct evidence of the formation of a 
transient accretion disk near the periastron passage. In 
the beat frequency model for QPO the Keplerian orbital 
frequency of the material at the inner edge of the accre- 
tion disk should exceed the spin frequency of the neutron 
star. Cheng & Ruderman (1989) have shown that in such a 
circumstance an electrostatic gap is open in the magneto- 
sphere around the null surface determined by f2» • B = 0. 
This gap can accelerate protons up to energies of tens of 
TeV, producing a hadronic current that impacts into the 
accretion disk generating tt° gamma-rays. The transient 
character of the disk makes the high-energy gamma radia- 
tion highly variable, as was found in the observed gamma- 
ray flux evolution. A specific prediction of the model is the 
suppression of gamma-ray emission when the column den- 
sity of the disk exceeds a critical value, near the peak of 
the X-ray luminosity. Future GeV and TeV observations of 
this source with instruments of high temporal resolution, 
like GLAST or 5@5 (see Aharonian et al. 2001), could be 
used to test the proposed model and, if it is basically cor- 
rect, to probe the evolution of the matter content on the 
accretion disk in this extraordinary X-ray binary. 
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